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An effective description of a general class of stochastic phase oscillators is presented. For this, the 
effective phase velocity is defined either by invariant probability density or via first passage times. 
While the first approach exhibits correct frequency and distribution density, the second one yields 
proper phase resetting curves. Their discrepancy is most pronounced for noise-induced oscillations 
and is related to non-monotonicity of the phase fluctuations. 



GENERAL INTRODUCTION 

The phase is a theoretical concept lying at the heart 
of description of oscillatory dynamics. In a physical con- 
text a periodic phase variable provides a simplified de- 
scription of self-sustained oscillations shown by natural, 
synthetic or mathematical systems in their state space 
[EH. The phase description contains many character- 
izing physical properties associated to oscillatory motion 
as for example frequency and regularity of oscillations. 
Most importantly, smooth or impulsive coupling of inter- 
acting oscillators may be formulated in terms of phase 
dynamics 0, 0, H . 

Irregular features, interpreted as noise, may be present 
in oscillations. In many situations noise can be treated 
as a perturbation to oscillatory dynamics. In this case 
one can start with the phase description of noiseless os- 
cillations, and consider noise as a perturbation. How- 
ever, noise may be substantial in the sense that it in- 
duces oscillations in the system which would equilibrate 
otherwise. Such noise-induced oscillations may be quite 
coherent disguising the underlying systems excitable na- 
ture. Therefore, it may not be possible to distinguish 
between the two cases in an experimental setup unless 
noise may be eliminated. 

Despite of evident similarities between noise-perturbed 
and noise-induced oscillations, a phase description of 
noise-induced oscillations cannot be obtained through 
standard perturbative procedures, because in the noise- 
free situation there is no dynamics. Recently, this prob- 
lem was addressed by an effective description, where 
noise-induced oscillations driven by an external periodic 
force could be characterized in a genuine way [BJ. It was 
seen, that the effective phase model relies on an average 
concept of speed given by the current velocity . 

The goal of this paper is to generalize the method 
of For this, we show that an alternative effective 
phase model is possible, and compare the two approaches. 
After a summary of well-known dynamical properties of 
stochastic phase oscillators in the next section, we outline 
the current model of effective phase theory for a single 
oscillator (Section J. Hereon, an alternative model based 
on first passage times is presented in Section „ modeling 
other aspects of stochastic dynamics. In Section,, mathe- 
matical aspects of the theory in its limiting cases for small 



and large noise are discussed highlighting differences and 
similarities of the effective phase models. A special em- 
phasis is put on the description of noise-induced oscilla- 
tions in these limits. In Section,, phase resetting curves of 
stochastic phase oscillators and their relation to the first 
passage model of effective phase theory are described. 

STOCHASTIC PHASE OSCILLATORS AND 
THEIR EFFECTIVE DESCRIPTION 

Our basic model is a stochastic phase oscillator. It 
is described by a 27r-periodic random process 9, called 
protophase, that obeys the Langevin dynamics 

9(t) = h(e(t))+ 9 (em(t), (i) 

where £(t) is ^-correlated Gaussian noise = 
25 (t — t'). A well-known example showing most promi- 
nent features of stochastic oscillations is the theta 
model 

6{t) = a + cos6 + a£{t) . (2) 

For \a\ < 1, it shows noise- induced oscillations which 
remain also in the deterministic limit a 2 — > 0, while for 
\a\ > 1 oscillations are noise-perturbed and they persist 
also for vanishing noise. 

For the stochastic phase oscillator introduced above, 
most relevant quantities are accessible analytically. The 
probability density P(6) is governed by the Fokker- 
Planck equation associated to equation (TTJ) given by 

d t P = -de [hP] + do [gdg [gP]] = -d e J . (3) 

For stationary probability density, the probability flux J 
is constant and we obtain the simpler equation 

J = hP-gd e [gP] . (4) 

This equation has the well-known solution Q 

PM-c/^-jS* (5) 
Je 0(%W>) 

where C is a normalization constant ensuring 
Jg 71 P(0) ol6 = 1. As illustrated in Fig. [TJ proba- 
bility density becomes singular for vanishing noise if 
oscillations are noise-induced. 



FIG. 1. Left: Probability density (JSJ of the theta-model 
for a — 0.9 and values of a as indicated. Right: Lyapunov 
exponent (|9]) of theta-model as a function of a for a — 0.5 
(blue circles), a — 0.95 (red triangles), and a — 1.5 (black 
squares) . 



FIG. 2. Mean frequency (|7|) in left plot, and quotient of 
diffusion coefficient ([8]) and mean frequency in right plot of the 
S-model as functions of a for a — 0.5 (blue circles), a = 0.95 
(red triangles), and a = 1.5 (black squares). 



Among the traditional quantities of stochastic phase 
oscillators are the mean frequency and diffusion coeffi- 
cient 
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which actually are properties of the point process n(t) 
counting the number of rotations of 9(t). They are ex- 
pressed in terms of h{9) and g{9) through the well-known 
formulas [H, Q 
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For noise-induced 



oscillations the mean frequency converges to zero in the 
limit of vanishing noise, as shown in the left plot of Fig. [3] 
In the large noise limit a 2 — > oo the mean frequency 
converges to a finite value for additive noise, which will 
be shown in a general setup in section „ The quotient of 
diffusion coefficient and mean frequency is a measure for 
decoherence of oscillations (right plot). For noise- induced 
oscillations, the well-known effect of coherence resonance 
is observable where decoherence decreases for increasing 
noise amplitude (red triangles and blue circles) [Io[ . The 
effect is not observable for noise-perturbed oscillations 
(black squares). 

Another interesting quantity that characterizes phase 
oscillations is Lyapunov exponent A associated to noise. 
It quantifies whether oscillators under the influence of 
the same noise representation £(t) will synchronize by 
stochasticity. For oscillator ([lj it is computed by 



\ = (h'(9)+g"(O)g(0)). 



(9) 



For oscillators under the influence by additive noise like 
the theta model, A vanishes in the limit a 2 — >• oo. For 
small noise there are three cases as exemplified in fig- 
ure ©. For a noise-perturbed oscillator (black squares) 
the Lyapunov exponent goes to zero as A cx —T For 
excitable oscillators (red triangles and blue circles) it con- 
verges to a finite negative value dependent on the quanti- 
tative stability of their fixed point. For these, there arises 
yet the special case where A shows a local minimum as 
a function of a (red triangles) which somehow resembles 
the noise-perturbed case. 

As we have seen, phase dynamics of noise-induced os- 
cillators can be quite close to that of periodic ones, what 
suggests a purely deterministic phase equation of the type 



= F{9) 



(10) 



for a theoretical description to be possible. Generally, 
one cannot set F{6) — h(0) by taking the deterministic 
part of the Langevin model ([1]), as h{0) may have zeros 
and thus 9 would be non-oscillating. Instead, we have 
to construct an effective phase model using some criteria 
to determine F, Generally, we demand that the effective 
phase model of type (fTU)) represent as many characteris- 
tic properties of stochastic phase oscillators as possible. 
In [5[ , one approach was proposed that allowed for a con- 
struction of an effective phase model with the same mean 
frequency (i) (equivalently, the same mean period) and 
distribution density (ii) as the stochastic phase oscillator. 
Without violating these conditions, the diffusion coeffi- 
cient (iii) could be modeled, too, by adding noise to (fTU)) 
in a certain way. Drawing a more general framework, we 
will present in the next two sections the effective phase 
model from Q and another effective phase model based 
on first passage times in a coherent way. Roughly, the 
difference of these models is in the definitions of the mean 
velocity. Given an interval of length A9, velocity can be 
measured as the quotient of A9 and the mean time At, 
that 9{t) spends in this interval. This leads exactly to 
an effective phase description as presented in [H}. Alter- 
natively, velocity can be measured as the quotient of A9 
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and the mean time taken to reach the opposite boundary 
of the interval, leading to the concept of first passage ve- 
locity. For deterministic phase oscillators the two defini- 
tions of velocity coincide, whereas for oscillations driven 
by a random force there is a difference, which is especially 
pronounced for oscillations that are noise-induced. 



CURRENT MODEL OF EFFECTIVE PHASE 
DYNAMICS 

We start by constructing the deterministic phase equa- 
tion for the current model 



9 = H(9) , 



(11) 



by using the notion of speed based on the mean time that 
the stochastic phase oscillator (flj spends in an interval 
[0,9 + d9]. This time is defined according to the invari- 
ant probability density dt — Pd9/J. It follows, that the 
current velocity H is given by 



H(8) 



Jd9 
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This leads to the same effective model introduced in [5|. 
Model ((TTj) obeys conditions (i) [preservation of the mean 
frequency] and (ii) [preservation of the probability den- 
sity] because it fulfills the stationary Liouville equation 
for P{9) and it shows correct period T as seen by 



d e [H{9)P{9)\ =0, and 
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In Fig. [3] we compare time series of oscillator ^ and its 
current model (fTTj) . Although the stochastic components 
are missing in the effective model, the observed dynamics 
is comparable. 

The current velocity H(9) can be expressed in terms 
of h{9) and g{9). For this, equation @ is divided by P, 
and the result is compared to equation (1121) yielding 



H(9) = h(9)--[g 2 {9)}' -g 2 {9) [lnP(9)]' = h{9)-u{9) . 

(13) 

It consists of the deterministic contribution h(9) and 
an osmotic contribution u(9), that is especially pro- 
nounced for the noise-induced oscillations in the theta 
model (Fig. [3] bottom panel). The current velocity 
corresponds to the point-wise average of central differ- 
ences 12]. Therefore, it may be constructed from an 
observed (e.g., experimentally) time series 9 n = 9(nAt) 
by a simple averaging procedure 



H{9) 



2At 



(14) 



Having constructed the current model, we can trans- 
form the protophase 9 to a uniformly rotating phase vari- 
able (/?, that describes oscillations in an invariant way. It 
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FIG. 3. Current model (|lip has dynamics 9(t) strongly re- 
sembling that of the theta model ([2}, here at a = 0.9 and 
a — 0.3 (top). For noise-induced oscillations the osmotic con- 
tribution u{9) to current velocity (|13p is non-negligible and 
accounts for its strict positivity. 



has simple properties 



f = uj; 
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As it can be easily checked, the nonlinear transformation 
f? — > is given by 



tp = S(9) = 2tt / P{9') d9' 
Jo 



(16) 



With the transformation, coordinate dependent differ- 
ences in the protophase A9 — 9 2 — 9\ can be transformed 
to invariant differences in the phase given by 

Aip = S{9 2 ) - S(0i) = 2tt f 2 P{9) d9 . (17) 

Given a set of data 9 n = 9{nAt) containing TV data 
points, transformation can be obtained numerically as 
in [3j. If one is not interested in the transformation S(9) 
but in the transformed data tp n — S(9 n ) only, one may 
alternatively evaluate 



2n N ~ 1 
(=0 



0, 



(18) 



which is implemented quickly by sorting (0 is the Heav- 
iside function). 

It was also shown in 0], that simultaneously to the 
mean period and the probability density, phase diffu- 
sion may be modeled. For this, (^-correlated noise rj{t) 
is added to the invariant phase dynamics (|15p with noise 
intensity \f~D leading to 



tp - 



Dr)(t) 
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Now, tp has a diffusion constant D while preserving uni- 
form density and mean frequency uj. Therefore, applica- 
tion of the inverse transformation — 5~ 1 (</?) gives us 
the stochastic current model 



= H(9) 



H(9) V (t) 



(20) 



that fulfills conditions (i) and (ii) for any value of D. 
It may be chosen freely, and we chose it uniquely from 
condition (iii): Diffusion coefficients of stochastic current 
model (|2"0"|) and stochastic phase oscillator (flj , should be 
equal. This condition is fulfilled if diffusion coefficient 
© is used for D. 



PASSAGE TIME MODEL OF EFFECTIVE 
PHASE DYNAMICS 

As an alternative to the current model outlined in last 
section, a velocity based on first passage time statistics 
of oscillator (JT]) shall lead us to a first passage model 



= N{6) 



(21) 



To determine 



A9 
At 



we interpret At as the first passage 



time of passing interval AO. More precisely, the first pas- 
sage velocity is constructed using the mean first passage 
time T(a, (3) which it takes for 9{t) to reach a boundary 
j3 > a starting at a. We get 



1 dt T(9,t 
— — - = — = hm 

N(9) d9 s^o 
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It fulfills condition that the mean frequency of (|2lj) 
should be equal to the mean frequency of the oscillations, 
because the mean period T is nothing else as T(0, 2ir + 9) 
and therefore f * d9/N{9) — T. However, the condition 
(ii) above [density in the effective model is equal to the 
density of stochastic oscillations] is generally not fulfilled, 
as will be shown next. 

The distribution density of model (|21[) . which we will 
call first passage density, is given by 



R{9) 
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In order to derive an equation for R(9), an equation for 
the mean first passage time has to be established Q. 
Consider Fokker-Planck equation Q, with the sharp ini- 
tial condition P(9,0) = 8(9 — a). In this case equation 
© describes the conditional probability P(9, t\a, 0). The 
boundary conditions 



P{-oo,t\a,0) = P(/3,t\a,0) = , 



(24) 



are introduced, which correspond to the fact that trajec- 
tories starting at 9(0) — a should only be considered as 
long as they do not reach boundary (3. Now, P has to 



be reinterpreted since the normalization condition does 
not hold anymore. The no-passage probability G(a, t) is 
defined as the probability that at time t boundary /? is 
not reached when starting at a. It is defined as 



G(a,t) 



P(9',t\a,0)d9' . 



(25) 



By a backward-Kolmogorov expansion of the Markovian 
transition probability P it is derived [|| that G obeys 

d t G = h(a)d a G + g(a)d a [g(a)d a G] . (26) 

For t £ [0, oo) and a < /3, a distribution density for 
random first passage times is given by g(a,t) = —dtG. 
With respect to it, the mean first passage time is given 
by 



T(a,P) = {t)= / G(a,t') dt 



(27) 



Integrating equation (|26p over positive times one obtains 
an equation for the mean first passage time 



-l = hd a T + gd a [gd a T] 



(28) 



Because d p T(9,9) = -<9 Q T(0,0) = l/N(9), equation 
(f28| may be rewritten for the first passage density R(9) 
as 



J = h,R + gdg [gR] , 



solved by 



R(9) = C 



dip 
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dip 



(29) 



(30) 



with a normalization constant C, and J — ui/2ir. Note 
that, equation (|2"9"|) is similar to but has an opposite 
sign of the second term. It provides an easy to handle 
analytic formula for the first passage velocity N(9). 

The passage time density has a direct meaning for the 
stochastic protophase 9(t), as we will explain in the fol- 
lowing. Let t n be the times of first passage, for which 
9(t < t n ) < 9(t n ) holds. This is illustrated in left 
panel of Fig. @] where realization 9(t) and point process 
9 n = 9(t n ) are counterposed. Because of the Markov 
property, 9 n gives the starting point for a measurement 
of passage time ending when 0(f) reaches 6 n +i = 9 n + d9. 
Although the trajectory of corresponding time segment 
0(t„ < t < t„+i) lies in the whole region < 0„ + i, for 
the first passage density R(9) it is attributed to the in- 
terval 9 n < 9 < 9 n+ i. So, in fact R(9) is the density of 
the envelope of the stochastic protophase 0(f), shown as 
a red bold line in Fig. |4] As shown in the right panel, 
the probability density P(9) and the first passage density 
R(9) can be quite different. The fact that R(9) is prob- 
ability density of the envelope of 0(f) can be used for a 
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numerical reconstruction from data, as it is shown later 
in this section. 

We would like to mention here that going from the 
stochastic protophase 6{t) to its envelope, we achieve 
a monotonically growing protophase. Indeed, phase is 
often understood as a strictly monotonic variable, in 
some sense a "replacement" for a time variable. For 
a stochastic oscillator one often observes "reverse" 
variations. Thus, taking the envelope is a natural way 
to restore a monotonic function of time. 




50 100 150 0.2 0.4 0.6 



t P{6),R{9) 

FIG. 4. Protophase 6(i) (black line) of theta model J2| and its 
envelope (red thick line) have similar long term dynamics (left 
plot), but their distribution densities P{0) (black circles) and 
R(0) (red squares) are different, especially for noise-induced 
oscillations, shown at parameters a — 0.9 and a = 0.3. 

The first passage model (|2"Tj) provides effective phase 
dynamics alternative to the current model. It fulfills 
condition (i) in that it shows the same mean frequency 
as oscillator ([1]), but instead of modeling probability 
density (ii), preserves first passage density (iib), which is 
for deterministic oscillators equal to distribution density. 
In Fig. [5j first passage (blue circles) and current velocity 
(red squares) are counterposed for theta models jl}, for 
cases of noise-induced (top panel) and noise-perturbed 
oscillations (bottom panel). In the latter N(0) and H(6) 
coincide for vanishing noise, whereas for noise-induced 
oscillations the difference widens. 



In the numerical example, current and first passage 
velocities are mapped to each other by mirror symme- 
try. This is due to the fact that both h and g are sym- 
metric in the theta model. For a stochastic phase oscil- 
lator with with symmetric h and g, the transformation 
9^—0 transforms equation ((29)) in (j4]). Their solutions 
are mapped to each other too. Therefore, symmetry of h 
and g implies N(6) = H(—6), as observed in Fig. [5] 

For the first passage model, a uniformly rotating phase 
variable is constructed by the transformation 

o 

4' = Z{6) = 2n J R{6') dO' . (31) 
o 
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FIG. 5. Current velocity (|12l) (red squares) and first passage 
velocity (1221) (blue circles) for the theta model ((2} differ es- 
sentially, for both noise-induced at a = 0.9 (top plot), and 
noise-perturbed at a = 1.1 (bottom plot) oscillations. Noise 
amplitude is a — 0.2. The difference in effective velocities is 
especially pronounced for noise-induced oscillations and does 
not disappear for small a 2 . 

Again, differences in protophase 82 — 9\ can are trans- 
formed to differences in phase by 

= Z(9 2 ) - Z{0i) = 2tt [ R(6) d6 . (32) 

As for the current model, noise may then be taken into 
account by 

= N{6) + —N(9)r](t) . (33) 

This stochastic first passage model fulfills conditions (i), 
(iib), and (iii). However, it does not fulfill condition (ii), 
as its probability density is R(0) 7^ P{6)- Therefore, 
an application of transformation (|31[) to the stochastic 
variable 9(t) does not lead to a uniformly distributed 
phase ip(t) as in equation (fT5j). 

The fact, that R(0) is the probability density of the en- 
velope of a realization 9(t), can be used to obtain trans- 
formation Z{ff) from a dataset 0(nAt). The envelope of 
data has to be constructed. For this, we need to find 
the times of first passages t n for which data has a his- 
tory that is strictly smaller, i. e. 9(jAt < t n ) < 8(t n ). 
At these first passages 6 n = 9(t n ), transformation (|3Tj) is 
estimated by formula 

To find the transformation on the whole domain one can 
proceed by an appropriate interpolation of Z(9 n ) as for 
example smoothing splines. Note, that equation (|3~4")l 
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gives a biased estimator. This can be fixed by insert- 



ing a central difference scheme (tj 



U 



0/2. 



ASYMPTOTIC PROPERTIES 
Singular perturbation for small noise 

In the Fokker-Planck equation ([3]), noise amplitude 
g(9) appears in front of the derivative with highest or- 
der in 9. Therefore, a perturbation expansion in g{9) for 
a deterministic approximation of dynamics ([TJ might be 
singular. For an oscillator with additive noise 



KB) + a^t) , 



(35) 



we show that perturbation expansion becomes singular 
if oscillations are noise-induced. We discuss the case in 
which the approximation of dynamics (UJ should model 
the distribution density of oscillator (f3"5j) . The singularity 
that arises in the limit a — > is discussed by perturbation 
expansion starting from a = 0, and starting from a > 
by effective phase theory. For the latter, current model 
jTTJ is used. 

For a = 0, our considerations start with the arbitrary 
model 9 = h(9) for oscillator (|35|) . Suppose, it shows 
noise-perturbed oscillations, i. e. h(9) strictly positive. 
Then, the distribution density of arbitrary model gives 
a proper zeroth order approximation to the probability 
density 



P{9) = 



C 

W) 



0(a) 



Better approximations can be obtained analytically by 
Taylor expansion of P(9). On the other hand, oscillator 
(|35p may show noise-induced oscillations, for which we 
impose without loss of generality that h(9) shows two 
zero crossings. Then, the model has a stable fixed point 
at and an unstable at 9 + > and its distribution 
density is given by 



(36) 



P(0) =6(6-0-)+O(<t) 



corresponds to the probability density of equation (|35|) . 
However, higher order terms in a that should lead to a 
smooth distribution density must necessarily be singular. 
It is seen, that perturbation theory becomes singular for 
noise-induced oscillations. Note, that first passage den- 
sity also shows a singular limit for noise-induced oscilla- 
tions, however, not at the stable, but the unstable fixed 
point. In the above example it is given by 



R(9) = 6(9 - 9+) + 0(a) 



(37) 



On the other hand, first passage density and probability 
density converge for noise-perturbed oscillations, reflect- 
ing the fact that for deterministic self-sustained oscilla- 
tors the quantities are synonym. 



Starting from a current model to oscillator (13S1 com- 
puted at finite <r, another view can be gained on the limit 
a —} 0. With formula (0 and (fl~2j) . current velocity cor- 
responding to ([35"1) is given by 



H{9) 
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+ 2TT dip 
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Je 

(38) 

Let us assume without loss of generality that r(0, 27r) 
is positive. In the limit of small noise, the integral in 
the denominator of equation (J38J) is dominated by the 
minimum of r(9, tp) with respect to ip. Using the method 
of stationary phase the zeroth order expansion 



H(9)=0(a) 



f H0) 
1° 



if argmin^ [r(9, ip)} = 9 



elsewise 



(39) 



is obtained. It can be seen, that the deterministic limit 
cannot be used at a = 0. Surprisingly, small noise limit 
leads to a vanishing of current velocity in a finite interval 
around the fixed point (9_ and not just at this point. A 
similar derivation is possible for first passage velocity, as 
for example seen by mirror symmetry. 

In Fig. [6] the peculiar nature of H(9) is illustrated for 
the theta model. If noise is small, the effective veloc- 
ity becomes discontinuous for noise-induced oscillations 
where \a\ < I (right plot), whereas for noise-perturbed 
oscillations at \a\ > I, H(9) converges to h(9) for all 9 
(left plot). 
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FIG. 6. Effective velocity at small a (red squares) converges 
to h (dashed black line) for noise-perturbed oscillations (left 
plot), but for noise-induced oscillations becomes discontinu- 
ous (right plot) as in formula (139|l . Furthermore, the total 
area under the curve is independent of a (formula (|42|l h Red 
squares: a = 0.1, Blue circles: a = 0.7. Parameter a is chosen 
to be 1.1 (left plot), and 0.9 (right plot). 



Estimating frequency for large noise 

For large noise, an estimation of frequency can be trou- 
blesome when one has to rely on Monte-Carlo simula- 
tions. Here, we want to provide a formula that allows 
for an estimation of frequency at large noise amplitude 
a — Y co , when given an effective velocity at arbitrary a. 
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In the limit of large noise the probability density be- 
comes uniform. By the definition of, for example, current 
velocity ([12"]) it is seen that: 



lim I 



lim — 

cr— >oo 2 



i r 2 * 

- / H(9) dd = lim 

In Jq er— S-oc 



(40) 



Now, it is shown that integral / does not depend on a. 
Using equation (fl"3)) . current velocity is represented as 
H = h — a 2 [In P]' which yields for the integral d a il = 0. 
Therefore, / does not depend on a, and we have 



i r 2n 

— / H{9) d9 = u 
W Jo 



(41) 



where H may be given for arbitrary noise intensity a. 
With equation (|39p it can also be verified that 



1 f 27T 
— / h{9) dd = L0 C 
^ Jo 



(42) 



For additive noise it is seen that the area under the func- 
tion H(9) is a preserved quantity and equal to 27ru; 00 . 
The same result can be obtained for the first passage 
velocity, straight-forwardly. 



STOCHASTIC PHASE RESETTING 

In the realm of phase resetting one is interested in the 
response of an oscillator to a brief stimulus, a kick, ap- 
plied to the system at certain protophase a with certain 
strength and direction k. The freely rotating period T is 
compared to the period Tk(a) wherein the kick is applied, 
by computing the phase resetting curve [l3| 



Ai/>(a,k) = 2?r 



T - T k {a) 
T 



(43) 



It gives the shift in uniformly rotating phase the oscil- 
lator experiences due to the kick. For example equation 
(|17p provides the phase resetting curve of current model 
([TTj) , whereas (|32|) provides it for first passage model (|2"Tj) 
when for each setting 62 = 9\ + k. Up to here, quantities 
appearing in equation (1431) are only well-defined for de- 
terministic systems that show limit-cycle oscillations. In 
this section we extend the applicability of formula (|4"3"[) 
to stochastic phase oscillators (cf. [3]) and discuss the 
results in terms of effective phase theory. 

A kick applied at a time t' with scalar strength k rep- 
resenting a brief stimulus is introduced to our stochastic 
phase oscillator (jTJ) by formula 

9(t) = h(9(t)) + g(9(t))t(t) + kS(t' - t) , (44) 

Phase resetting is computed by a comparison of the 
kicked and unkicked stochastic phase oscillator. For this, 
quantities appearing in formula (|43[) are interpreted as 



follows: Quantity T is given by mean period T = 2n/to 
computed for the unkicked oscillator. Quantity Tk(a) is 
computed for the kicked oscillator as the mean condi- 
tional first passage time starting at 9 = a + k (value just 
after kick), to reach boundary 9 = a + 2tt. Here, the 
mean shall be taken with respect to noise. 

By the above interpretation of formula ([43]) for kicked 
stochastic oscillators, phase resetting can be calculated 
using mean first passage times. From formula (|28|) it is 
deduced that T(a, /3) + T(f3, 5) = T(a, 5). This is related 
to the Markov property of 9(t), and it allows us to express 
mean conditional first passage time as 



Tk(a) = T — T{a, a + k) , 



(45) 



For equation (j43]), it follows Aip(a, k) = 2nT(a, a+k)/T. 
This can be rewritten in terms of R using equation 
as 



AiP(a,k) = 2tt J R{9) d9 



(46) 



This is the exact formula of phase resetting curve 
Atp(a, k) for a general stochastic phase oscillator (jTJ) . 
It is seen that, phase resetting curve 



a+k 



A 



V (a, k) = lim [<ptf + E ) - tp{H - e)] m=a = 2n J P(9) 

a 

(47) 

derived for the current model does not correspond to that 
of the stochastic phase oscillator, whereas phase resetting 
curve 



d9 



a-\-k 



AV(a, k) = lim W + e) - W - e) V)= a 



= 2tt 



P{9) d9 



(48) 

derived from first passage model ([21]) yields the correct 
formula (|46j) . Let us explain why the current model 
fails. In section „ it was seen that current velocity is 
given by H(9) — J/P, which leads to the stationary 
solution (fT2| . However, the stationary state is broken 
in the phase resetting procedure, where the time evo- 
lution of 9(t) starts from the definite value 9 = a + k. 
Consider for example the moment, right after the reset- 
ting, where P(9) = 6(6 — a — k). The probability flux 
is calculated by integrating equation and one obtains 
J{a) = h(a + k). Deviations to the stationary solution 
(|12|) are most prominent in the excitable regime where 
J(a+k) = h(a+k) < 0. In this case, the time-dependent 
current model 



J(M) 
P(9,t) 



(49) 



has non-monotonic dynamics, and therefore does not 
yield a good phase description. 
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The failure of current model to predict the correct 
phase resetting curve is illustrated in figure [7] by noise- 
induced oscillations of the theta model. For this, 1000 
representations starting from 9% = 9 + — 0.2 (red solid 
line) and 62 = + + 0.2 (blue dashed line) were cal- 
culated and transformed to uniformly distributed phase 
(p = S(6). As before, 9+ is the unstable fixed point of 
corresponding arbitrary model. For each group an aver- 
age over representations was performed to obtain average 
dynamics tp(t). Equation (|15|) predicts tp to be rotating 
uniformly. Numerically, it is found however that while 
asymptotically uniformly rotating, initially tp(t) has sys- 
tematic non-uniformities such that initial phase differ- 
ence Aip(0) — S(92) — S(9i) widens. Asymptotically, it 
reaches the value Aijj = Z(02) — Z{9{) as correctly pre- 
dicted by stochastic phase resetting and first passage 
model (|48|) . For the stochastic phase oscillator with ini- 
tial condition 9\ with h(9\) < 0, it is furthermore seen 
that 37 (y) < initially, as explained above. Performing 
the same procedure using transformation ip = Z(9) it is 
confirmed that formula (|48p predicts the correct phase 
resetting curve (see right plot). 



CONCLUSIONS 

Two useful formulations of effective phase theory are 
presented, each relying on different definitions of phase 
velocity. While the concept of time spent on average in 
an interval leads to current velocity, a definition based 
on first passage times leads to first passage velocity. For 
noise-perturbed oscillations, the two velocities converge 
for vanishing noise, whereas for noise-induced oscillations 
they don't. Whereas the current model was shown to be 
useful for a characterization of continuous coupling [5j 
of stochastic phase oscillators, it was seen in this article, 
that first passage model gives a correct description of 
phase resetting. 

A natural extension of effective phase theory would 
be to construct models of two-dimensional oscillators 
for which a phase variable and a corresponding radius 
variable may be found. The extension is not straight- 
forward: A deterministic current-type model for higher- 
dimensional oscillators must have a meaning different 
from the one-dimensional case, because the correspon- 
dence between distribution and probability density (con- 
dition (ii)) can not be drawn. Here, also other interest- 
ing difficulties arise as described in [ljj . On the other 
hand, it might be possible to construct a first passage- 
type model exhibiting the correct phase resetting curve. 
In this sense the approach already implies, that a pertur- 
bative approach, as applied in [16j . for a construction of 
a phase resetting curve is expected to be biased. 

Together with these issues, it remains unresolved how 
to construct a theoretical description of a stochastic 
phase oscillator under the influence of an external force 





FIG. 7. Dynamics of tp(t) = S(0(t)) and tp(t) = Z(6(t)), 
averaged over 1000 representations starting at two distinct 
initial conditions B\/2 = 9+ ± 0.2 (top plots, +: blue dashed, 
— : red solid). Left column confirms that ()15[) does not hold 
for the phase resetting. For fixed initial conditions the initial 
phase difference A<^(0) = S(02) — S(8i) widens (bottom plot, 
top Fig.) to the value Alp = Z{t&2) — Z{9\) (red horizontal 
line); for larger times phase difference diffuses away due to 
insufficient averaging. Right column confirms that equation 
(|48[) gives the correct phase resetting curve for stochastic os- 
cillators. The initial phase difference Aip(0) = ZiQ^) — Z(9\) 
persists (bottom plot, bottom Fig.) at its theoretically pre- 
dicted value. 



consisting of both, a continuous and a pulse-like contri- 
bution. Furthermore, we expect the presented work to 
be useful for a characterization of finite ensembles of os- 
cillators with either continuous or pulse-like coupling. 
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